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This  paper  describes  methods  and  results  of  research  for  incorporating  four  different  parameterized 
wave  breaking  and  dissipation  formulas  in  a  coastal  wave  prediction  model.  Two  formulations  assume 
the  breaking  energy  dissipation  to  be  limited  by  the  Rayleigh  distribution,  whereas  the  other  two 
represent  the  breaking  wave  energy  by  a  bore  model.  These  four  formulations  have  been  implemented 
in  WABED,  a  directional  spectral  wave  model  based  on  the  wave  action  balance  equation  with 
diffraction,  reflection,  and  wave-current  interaction  capabilities.  Four  parameterized  wave  breaking 
formulations  are  evaluated  in  the  present  study  using  two  high-quality  laboratory  data  sets.  The  first 
data  set  is  from  a  wave  transformation  experiment  at  an  idealized  inlet  entrance,  representing  four 
incident  irregular  waves  in  a  slack  tide  and  two  steady-state  ebb  current  conditions.  The  second  data  set 
is  from  a  laboratory  study  of  wave  propagation  over  a  complex  bathymetry  with  strong  wave-induced 
currents.  Numerical  simulation  results  show  that  with  a  proper  breaking  formulation  the  wave  model 
can  reproduce  laboratory  data  for  waves  propagating  over  idealized  or  complicated  bathymetries  with 
ambient  currents.  The  extended  Goda  wave  breaking  formulation  with  a  truncated  Rayleigh 
distribution,  and  the  Battjes  and  Janssen  formulation  with  a  bore  model  produced  the  best  agreement 
between  model  and  data. 

©  2008  Published  by  Elsevier  Ltd. 


1.  Introduction 

As  waves  approach  the  shore,  their  heights,  lengths  and 
directions  are  changed  due  to  shoaling,  refraction,  diffraction, 
reflection  and  breaking  as  a  consequence  of  the  particular 
bathymetric  and  geometric  features,  as  well  as  encounter  with 
currents  and  structures  present.  Because  reliable  wave  predictions 
in  the  coastal  areas  are  crucial  to  engineering  applications 
associated  with  shore  protection,  sediment  management,  harbor 
construction,  navigation  channel  maintenance  and  maritime 
disaster  reduction,  numerical  modeling  of  nearshore  wave 
transformation  has  been  a  subject  of  considerable  interest.  In 
the  past  three  decades,  significant  advances  have  occurred  in  the 
wave  modeling.  Many  wave  transformation  models  have  been 
developed  by  using  different  type  of  linear  and  nonlinear  wave 
theories  for  propagation  of  monochromatic  or  irregular  waves 
over  an  arbitrary  bathymetry.  Because  all  wave  models  have 
certain  advantages  and  limitations,  their  appropriateness  in  the 
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coastal  applications  is  largely  dependent  on  type  of  physical 
processes  at  the  project  site. 

Because  waves  are  random  in  the  real  world  setting,  it  is 
necessary  to  take  into  account  their  random  nature  in  the  wave 
transformation  models.  Mase  and  Kitano  (2000),  Nwogu  and 
Demirbilek  (2001),  and  Zubier  et  al.  (2003)  have  classified  random 
wave  transformation  models  into  two  categories.  The  first  category 
includes  models  based  on  the  energy  balance  equation  or  wave 
action  equation.  Examples  are  TOMAWAC  (Benoit  et  al.,  1996), 
GHOST  (Rivero  et  al.,  1997),  STWAVE  (Smith  et  al.,  1999),  SWAN 
(Booij  et  al.,  1999),  and  EBED  (Mase,  2001 ).  These  frequency-domain 
models  are  suited  to  directional  wave  transformation  over  large 
areas  of  open  oceans  but  have  in  recent  years  become  increasingly 
popular  for  modeling  neashore  waves.  The  second  category  of  wave 
models  includes  phase-resolving,  refined  time-domain  models 
based  on  the  conservation  of  mass  and  momentum  equations. 
The  Boussinesq  family  of  wave  models  belongs  to  this  class.  These 
phase-resolving,  time-domain  models  are  generally  more  resources 
demanding,  and  appropriate  for  applications  to  relatively  small 
coastal  areas  in  shallow  water  to  accurately  represent  wave  profile, 
wave  set-up,  wave-induced  current  by  solving  mass  and  momen¬ 
tum  equations  (e.g.  Peregrine,  1967;  Madsen  and  Sorensen,  1992; 
Nwogu,  1993;  Nwogu  and  Demirbilek,  2001). 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2QQg  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2008  to  00-00-2008 

4.  TITLE  AND  SUBTITLE 

Implementation  and  evaluation  of  alternative  wave  breaking  formulas  in 
a  coastal  spectral  wave  model 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

State  Key  Laboratory  of  Hydrology-Water  Resources  and  Hydraulic 
Engineering, Hohai  University, Nanjing,  Jiangsu  210098,  China, , 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

12 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


J.  Zheng  et  al.  /  Ocean  Engineering  35  (2008)  1090-1101 


1091 


In  some  coastal  areas  such  as  tidal  inlets  and  estuaries,  where  a 
great  deal  of  human  activity  is  located,  the  tidal  currents  can  be 
strong  and  their  effects  on  wave  transformation  cannot  be 
neglected.  Often  waves  are  shortened  and  steepened  by  ebb 
currents,  leading  to  considerable  wave  breaking  outside  and  in  the 
navigation  channels.  If  currents  are  strong,  wave  blocking  can 
happen,  and  considerable  navigation  hazard  conditions  can  occur. 
For  flood  tides,  waves  are  lengthened  and  their  heights  are 
reduced  by  current.  In  this  case,  larger  non-breaking  waves  inside 
the  channel  may  occur  as  compared  to  the  ebb  current  case.  In 
parts  of  the  nearshore  areas  with  a  complicated  bathymetry, 
strong  wave-induced  currents  can  create  converging  or  diverging 
wave  energy  zones.  The  Doppler  shift  phenomenon  affects  wave 
refraction,  reflection,  and  breaking  that  can  significantly  modify 
the  overall  redistribution  of  wave  energy  over  the  wave  spectra 
frequencies  and  directions.  Under  such  circumstances,  the  effects 
of  ambient  currents  on  nearshore  wave  transformations  must  be 
taken  into  account  in  the  wave  predictions  for  coastal  projects. 

One  of  important  factors  in  modeling  nearshore  waves  is  the 
choice  of  wave  breaking  formulation.  In  the  nearshore  zone,  once 
waves  start  to  break,  the  turbulent  dissipation  of  energy  becomes 
the  dominant  dissipative  mechanism,  and  breaking  processes 
control  the  spatial  variation  of  wave  heights.  In  the  absence  of 
ambient  currents,  Zhao  et  al.  (2001)  and  Zubier  et  al.  (2003)  have 
shown  that  different  formulas  for  parameterization  of  wave 
breaking  can  yield  large  differences  in  wave  height  estimates  in 
the  surf  zone.  They  demonstrated  with  a  two-dimensional  elliptic 
wave  model  and  a  spectral  wave  model  that  the  formulation  of 
Battjes  and  Janssen  (1978)  generally  produced  a  good  agreement 
between  numerical  estimates  and  data.  In  the  presence  of  the 
ambient  currents,  many  studies  have  been  conducted  to  investi¬ 
gate  effects  of  current  on  wave  breaking  (e.g.,  Yu,  1952;  Iwagaki 
et  al.,  1980;  Hedges  et  al.,  1985;  Lai  et  al.,  1989;  Sakai  et  al.,  1988; 
Li  and  Dong,  1993;  Briggs  and  Liu,  1993;  Ris  and  Holthuijsen,  1996; 
Smith  et  al.,  1998;  Chawla  and  Kirby,  2002).  However,  no 
comprehensive  evaluation  of  different  wave  breaking  formula¬ 
tions  in  two-dimensional  numerical  models  has  been  conducted 
for  random  wave  transformation  over  a  complicated  bathymetry 
with  ambient  currents.  A  field  evaluation  of  wave  breaking 
parameterizations  has  been  reported  by  Smith  (2001).  Lin  and 
Demirbilek  (2005)  used  a  data  set  for  irregular  waves  collected 
around  an  ideal  tidal  inlet  in  the  laboratory  as  a  benchmark  to 
examine  the  performance  of  two  spectral  wave  models.  They 
found  that  both  models  considerably  underestimated  wave 
heights  seaward  of  the  inlet,  suggesting  that  further  investigation 
of  wave  breaking  and  wave-current  interaction  near  inlets  are 
necessary  for  improved  spectral  wave  model  estimates  in  coastal 
inlet  projects. 

The  purpose  of  the  present  study  was  to  evaluate  four 
parameterized  wave  breaking  formulations  implemented  in 
WABED  (Mase  and  Kitano,  2000;  Mase,  2001;  Mase  et  al., 
2005a,  b;  Lin  et  al.,  2006,  2008)  for  coastal  spectral  wave 
transformation.  This  is  accomplished  by  adding  different  dissipa¬ 
tion  formulas  to  WABED  and  evaluating  estimate  of  wave 
dissipation  with  different  wave  breaking  formulations.  Model 
predictions  are  compared  to  data  collected  from  two  laboratory 
experiments.  In  the  first  experimental  study,  wave  measurements 
were  made  to  determine  effects  of  wave  shoaling,  breaking  and 
steady  ebb  currents  around  an  idealized  inlet  (Smith  et  al.,  1998). 
This  study  covered  a  wide  range  of  wave  and  current  parameters, 
and  data  obtained  can  be  used  in  evaluation  of  wave  dissipation 
formulations  for  current-  and  depth-induced  wave  breaking.  The 
second  experimental  study  investigated  random  wave  transfor¬ 
mation  accompanied  with  wave  breaking  over  a  complicated 
bathymetry,  where  strong  wave-induced  nearshore  currents  were 
occurring. 


2.  Model  description 

2 A.  Wave  action  balance  equation  with  diffraction  effect 


A  steady-state  spectral  model  WABED  (Mase  et  al.,  2005a,  b) 
is  used  to  evaluate  the  parameterized  wave  breaking  formula¬ 
tions.  It  is  a  2-D  phase-averaged  model  that  neglects  changes 
in  the  wave  phase  in  calculating  wave  and  other  nearshore 
processes  from  the  wave  energy  density.  Both  wave  diffraction 
and  reflection  are  included  in  approximate  ways.  To  take  the 
effect  of  ambient  currents  into  account,  the  wave  action  density 
is  used  in  WABED  rather  than  the  wave  energy  density  since 
the  wave  action  density  is  conserved  whereas  the  wave  energy 
density  is  not  if  waves  travel  with  ambient  currents  (Bretherton 
and  Garrett,  1968).  The  Doppler  shift  is  considered  in  the  wave 
dispersion  equation  implemented  in  WABED  to  develop  a 
practice-oriented  random  wave  model  for  coastal  engineering 
studies  at  inlets,  navigation  projects,  and  wave-structure  interac¬ 
tions  (Lin  et  al.,  2006,  2008).  In  these  applications,  wave  breaking, 
dissipation,  reflection,  diffraction,  and  wave-current  inter¬ 
action  are  always  important  processes,  and  accurate  representa¬ 
tion  of  these  processes  is  necessary  for  reliable  estimate  of 
waves  in  engineering  design,  maintenance,  and  operations.  This 
paper  investigates  effects  of  parameterized  wave  breaking  and 
dissipation  on  wave  transformation  in  inlets  and  navigation 
channels.  The  frequency-dependent  variation  of  wave  action 
density  resulting  from  wave-wave  interactions  is  not  considered 
in  the  present  study,  and  will  be  addressed  in  a  companion  future 
paper. 

The  governing  wave  action  balance  equation  with  the 
wave  diffraction  effect  as  implemented  in  WABED  model  (Mase, 
2001)  is 


6  (C*N)  d(CyN)  8(CeN) 
8  x  +  8y  +  80 
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2(7 


(CCgCOS2  0Ny)y  - 
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CCg  COS2  ONyy 


-ebN 
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where  N  is  the  wave  action  density,  defined  as  wave  energy 
density  divided  by  the  angular  frequency  a  relative  to  a  current 
(Doppler  shift).  The  horizontal  coordinates  are  (x,  y),  and  6  is  the 
wave  direction  measured  counterclockwise  from  the  x-axis.  The 
first  term  in  the  right-hand  side  of  Eq.  (1)  represents  wave 
diffraction  as  formulated  from  the  parabolic  wave  approximation 
assumption.  A  default  value  of  k  =  2.5  is  used  for  the  diffraction 
intensity  parameter  in  the  present  study.  As  suggested  by  Mase 
(2001),  the  appropriate  values  of  k  should  be  estimated  from 
laboratory  and  field  data,  and  in  the  absence  of  data,  the 
recommended  value  is  2.5.  In  Eq.  (1),  C  and  Cg  are  the  wave 
celerity  and  group  velocity,  respectively,  and  cb  is  a  parameterized 
wave  breaking  function  for  wave  energy  dissipation.  The  char¬ 
acteristic  wave  velocities  with  respect  to  x,  y  and  0  coordinates  are 
accordingly  Cx,  Cy,  and  C0,  defined  as 

Cx  =  Cg  cos  0  +  U  (2) 
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where  U  and  V  are  current  velocity  components  in  the  x  and  y 
directions,  respectively,  and  k  is  the  wave  number.  The  relation¬ 
ships  between  the  relative  angular  frequency  a,  absolute  angular 
frequency  at,  wave  number  vector  k,  current  velocity  vector  U,  and 


1092 


].  Zheng  et  al.  /  Ocean  Engineering  35  (2008)  1090-1101 


the  water  depth  h  are 

a2  =  g|k|  tanh  |k|h 

(5) 

n  —  cj  k  U 

(6) 

The  wave  action  balance  equation  with  diffraction  (Eq.  (1))  is 
solved  by  a  forward-marching  first-order  upwind  finite-difference 
method.  For  given  values  at  the  offshore  boundary,  wave  spectra 
and  statistical  quantities  are  calculated  in  the  wave  propagation 
direction  at  each  column  in  a  rectangular  grid  before  moving 
forward  to  the  next.  WABED  can  optionally  perform  the  backward 
marching  for  seaward  reflection  after  completing  the  forward¬ 
marching  calculations. 

In  the  treatment  of  the  dispersion  relation  with  the  Doppler 
shift,  waves  with  large  wave  numbers  were  not  considered  in  the 
original  WABED  (Mase  et  al.,  2005a).  However,  the  exclusion  of 
large  wave  numbers  (fc)  can  in  some  cases  lead  to  errors.  In  the 
present  version,  we  use  a  quadratic  solution  for  wave  dispersion 
parameter  kh>3i:,  and  set  tanhkh  =  1  for  large  kh,  and  combine 
Eqs.  (5)  and  (6)  into  one  quadratic  equation  in  terms  of  wave 
number  as 

k2(U  cos  0  +  V  sin  0 )2  -  [2 (U  cos  0  +  V  sin  0)co+g]k 

+  co2  =  0  (7) 

The  solution  for  the  wave  number  in  Eq.  (7)  is  the  smaller 
real  root  in  order  for  the  wave  group  velocity  to  be  greater  than 
the  ambient  current.  For  kh<3n,  the  solution  is  obtained  through 
an  iterative  scheme.  If  no  solution  exists  for  the  dispersion 
relationship  indicating  a  wave  blocking  condition,  the  wave  action 
density  is  set  to  zero  for  the  corresponding  frequency  and 
direction  bin. 

2.2.  Parameterization  of  wave  breaking  energy  dissipation 


et  al.  (1980)  verified  that  the  Miche’s  breaker  criterion  (1944) 
can  be  applied  to  a  wide  range  of  experimental  data  over  a 
sloping  beach  with  currents  if  the  wavelength  is  calculated 
with  current  included  in  the  wave  dispersion  relation. 
Battjes  (1972)  extended  the  Miche’s  criterion  for  variable 
depth  as 


^  =  0.14  tanh 
i-b 


(9) 


where  Lb  is  the  wavelength  at  the  breaking  point  with  currents 
considered,  and  yb  is  an  adjustable  coefficient  that  varies 
with  beach  slope.  In  application  to  irregular  waves,  a  constant 
value  of  yb  =  0.8  is  frequently  adopted  for  a  wide  range  of 
beach  slopes  and  mean  wave  steepnesses  (Battjes  and  Janssen, 
1978;  Zhao  et  al.,  2001;  Zubier  et  al.,  2003).  In  WABED,  a 
functional  form  as  suggested  by  Ostendorf  and  Madsen  (1979)  is 
adapted 


j  0.8  +  5  tan  /?  tan/?  <0.1 
|  1.3  tan  /?>0.1 


(10) 


where  /?  denotes  the  local  bottom  slope  angle. 

The  change  in  breaker  height  with  respect  to  the  cell  length  dl 
is  expressed  as 


d  Hb/dl  = 
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(J.OO 
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The  breaking  heights  at  each  side  of  a  cell  in  the  computational 
grid  are 


Hw  =  fib  —  2  dUb 


(12) 


In  Eq.  (1),  the  parameterized  function  cb  describes  the  mean 
energy  dissipation  rate  per  unit  horizontal  area  due  to  wave 
breaking.  The  importance  of  this  function  is  examined  for 
four  wave  breaking  formulas:  (1)  extended  Miche  formula, 
(2)  extended  Goda  formula,  (3)  Chawla  and  Kirby  formula,  and 
(4)  Battjes  and  Janssen  formula.  These  formulas  are  divided  into 
two  generic  categories.  The  first  category  describes  the  wave 
breaking  energy  dissipation  with  various  breaker  criteria  and  by 
truncating  the  tail  of  the  Rayleigh  distribution.  The  second 
category  describes  the  integral  (bulk)  wave  energy  dissipation 
over  the  spectrum  with  a  bore-type  model.  For  additional 
information,  see  Zhao  et  al.  (2001)  and  Zubier  et  al.  (2003).  A 
summary  description  of  the  parameterized  wave  breaking 
formulas  is  presented  next. 

2.2.1.  Extended  Miche’s  formula 

Following  Takayama  et  al.  (1991),  Mase  et  al.  (2005a) 
formulated  the  wave  energy  dissipation  rate  Eb  using  the  extended 
Miche’s  breaking  criteria: 


Hbo  —  Wb  + 


1 

2 


dHb 


(13) 


2.2.2.  Extended  Goda's  formula 

Based  on  the  laboratory  data,  Goda  (1970)  proposed  a  breaker 
criterion  by  taking  into  account  the  bottom  slope  and  deep  water 
wave  steepness.  In  the  presence  of  opposing  currents,  Sakai  et  al. 
(1988)  introduced  a  proportional  coefficient  to  account  for  the 
combined  effects  of  unit  width  discharge,  bottom  slope,  incident 
wave  length,  and  local  water  depth.  With  this  modification,  the 
breaking  wave  height  becomes 
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where  Hbi  and  Hbo  are  breaking  wave  heights  at  the  seaside 
and  landside  of  a  cell,  Hi/3  is  the  significant  wave  height, 
and  dl  is  the  grid  size.  Eq.  (8)  is  valid  if  Hbi>Hbo  and  eb  =  0 
otherwise.  The  extended  Miche’s  formula  reduces  to  a  steepness 
limit  in  deep  water  and  a  depth  limit  in  shallow  water.  Iwagaki 


(  0.506  Ed>  0.0024 

c(ed )  =  <  1.13  — 260ed  0.0024  >£d^  0.0005  (15) 

(  1.0  sd<  0.0005 

and 

Ed  =  ^3  (tan1/4/?)  / (7i/Lo)  (16) 

S  T  / 

where  L0  is  the  wavelength  in  deep  water,  q  the  unit  width 
discharge,  g  the  gravitational  acceleration,  and  T  the  significant 
wave  period  for  random  waves. 
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The  change  in  breaker  height  with  respect  to  the  cell  is  2.2.4.  Battjes  and  Janssen’s  formula 
defined  as  Battjes  and  Janssen  (1978)  developed  the  following  formula  for 

the  prediction  of  energy  dissipation  in  random  waves  breaking  on 


d  Hb/dl 


-0.225k  tan  /?( 1  +  15tan4/3  /J)c(£d)  x  exp^-1.5f^(l  +  15tan4/3/bj  tan  fi^O 
0  tan  /?<  0 


(17) 


The  energy  dissipation  rate  eb  is  calculated  by  combining  Eqs. 
(8),  (12)-(14)  and  (17). 


2.2.3.  Chawla  and  Kirby’s  formula 

Thornton  and  Guza  (1983)  parameterized  the  wave  breaking 
energy  dissipation  using  the  Rayleigh  distribution  and  an 
empirical  weighting  function  in  a  bore-type  model.  However, 
because  their  formula  is  independent  on  the  local  water  depth,  it 
cannot  be  used  to  model  current-limited  wave  breaking.  Using  a 
similar  approach,  Chawla  and  Kirby  (2002)  proposed  the  follow¬ 
ing  expression  for  bulk  energy  dissipation  with  the  ambient 
current: 


(D) 


3Xp  /W  /  k  \ 

32  y  tanh  kh  (y  tanh  kh J 


x  H 


5 

rms 


(  feHrms  \ 

l  y  tanh  kh  J 


-5/21 


(18) 


where  <D>  is  the  bulk  energy  dissipation  by  all  breaking  waves,  k 
is  the  wave  number  corresponding  to  the  mean  angular  frequency 
<7,  Hrms  is  the  root-mean-square  wave  height,  and  scaling 
parameters  l  and  y  are  set  to  0.4  and  0.6,  respectively. 

The  wave  breaking  energy  dissipation  coefficient  eb  in  Eq.  (1 )  is 
calculated  as 


(■ P ) 

(Pgrfm  s/8)^ 


(19) 


a  beach: 


(D)=afQ.jH't  (20) 

where  a  is  a  constant  of  order  one, /(=  a/2n)  is  the  mean  relative 
frequency,  and  Qt,  is  the  fraction  of  broken  waves  greater  than  the 
breaker  height  Hb  in  the  Rayleigh  distribution  and  as  determined 
from  the  following  equation: 


1-Qb_  (Hrms\2 
In  Qb  \Hb  J 


(21) 


In  SWAN  version  40.11  (Booij  et  al„  1999),  a  constant  breaker 
parameter  of  yb  =  0.73  was  used  in  Eq.  (9)  to  determine  the 
breaking  wave  height  Hb  at  the  local  water  depth.  Chen  et  al. 
(2005)  also  adopted  this  breaking  criterion  in  their  finite  element 
coastal  harbor  wave  model  based  on  the  extended  mild  slope 
equation  with  wave-current  interaction.  Eqs.  (20)  and  (21)  are 
used  with  Eq.  (19)  to  calculate  eb. 


3.  Numerical  results  and  discussion 

3.1.  Wave  shoaling  and  breaking  around  an  idealized  inlet 

Smith  et  al.  (1998)  conducted  a  laboratory  study  to  investigate 
the  wave-current  interaction  and  wave  breaking  in  the  presence 
of  a  steady  ebb  current  at  an  idealized  entrance.  The  physical 
model  was  built  at  an  undistorted  model-prototype  scale  of  1:50. 
The  shore-parallel  depth  contours  were  determined  by  an 
equilibrium  profile  equation  of  Dean  (1977)  using  a  coefficient 
of  0.24.  The  sloping  bathymetry  in  this  model  extended  to  18.2  cm 


Table  1 

Incident  wave  parameters  and  current  conditions  for  idealized  inlet  case 


Hm0  (cm) 

fp  (Hz) 

Have 

(cm/s) 

Notes 

Case  1 

5.77 

0.78 

11.5 

Run  5  in  the  Technical 

Case  2 

3.92 

1.45 

11.1 

Report  CHL-98-31 

Run  8  in  the  Technical 

Case  3 

5.97 

0.78 

21.9 

Report  CHL-98-31 

Run  9  in  the  Technical 

Case  4 

5.51 

1.37 

21.9 

Report  CHL-98-31 

Run  11  in  the  Technical 
Report  CHL-98-31 

a 

400 
200 
0 

0  200  400  600  800  1000  1200  1400  1600 


b 

400 
200 
0 

0  200  400  600  800  1000  1200  1400  1600 


moderate  opposing  currents 


strong  opposing  currents 


Fig.  2.  Input  current  fields  for  idealized  inlet  case:  (a)  moderate  opposing  currents  and  (b)  strong  opposing  currents. 
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0.05  n 


0.04  - 


0.03  - 


0.02  - 


0.01  - 


Hm(j=  5.77cm,  fp=0.78Hz,  U=  14cm/s 

•  Extended  Miche’s  formulation  (1979) 
Extended  Goda's  formulation  (1989) 
Chawla  and  Kirby's  formulation  (2002) 

•  Battjes  and  Janssen’s  formulation  (1 978) 


400  800  1200  1600  2000 

x(cm) 


Fig.  3.  Model  results  of  wave  height  and  energy  dissipation  along  central  transect  for  small  peak  frequency  waves  with  moderate  opposing  currents:  (a)  wave  height  and 
(b)  energy  dissipation  coefficient. 


Fig.  4.  Model  results  of  wave  height  and  energy  dissipation  along  central  transect  for  large  peak  frequency  waves  with  moderate  opposing  currents:  (a)  wave  height  and 
(b)  energy  dissipation  coefficient. 


Fig.  5.  Model  results  of  wave  height  and  energy  dissipation  along  central  transect  for  small  peak  frequency  waves  with  strong  opposing  currents:  (a)  wave  height  and 
(b)  energy  dissipation  coefficient. 
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Fig.  6.  Model  results  of  wave  height  and  energy  dissipation  along  central  transect  for  large  peak  frequency  waves  with  strong  opposing  currents:  (a)  wave  height  and 
(b)  energy  dissipation  coefficient. 


Fig.  7.  Model  to  data  comparisons  using  Chawla  and  Kirby’s  formulation  and  its 
modifications. 


Table  2 

Statistical  mean  errors  and  correlation  coefficients  for  idealized  inlet  case 


MAERH  (%) 

Correlation 

coefficient 

Extended  Miche’s  formulation 

18.9 

0.88 

Extended  Goda’s  formulation 

13.1 

0.95 

Chawla  and  Kirby’s  formulation  (y  = 

0.6) 

25.5 

0.79 

Chawla  and  Kirby’s  formulation  (y  = 

1.0) 

15.8 

0.91 

Chawla  and  Kirby  ‘s  formulation 

21.3 

0.85 

(y  =  0.6  add  breaker  criterion) 

Chawla  and  Kirby’s  formulation 

22.7 

0.84 

(y  =  0.6  add  breaker  percentage) 

Battjes  and  Janssen’s  formulation 

14.5 

0.93 

Section  x=9.0m 


Section 

y=6.0m 


Fig.  8.  Bathymetry  and  locations  of  wave  gauges  and  currentmeters  (circles) — 
depth  contours  shown  as  dash  lines  and  labeled  in  centimeters. 


mean  low  water  depth,  and  was  linearly  transitioned  to  the  basin 
floor  at  a  depth  of  30.4  cm.  Two  parallel  jetties  extending  550  cm 
offshore  were  spaced  366  cm  apart,  and  built  to  create  an  entrance 
channel  in  which  depths  varied  from  12.8  to  9.0  cm.  The  inlet 
throat  region  converged  to  a  depth  of  15.2  cm  relative  to  a  mean 
low  water  datum.  Incident  waves  and  current  conditions  were 
widely  varied  in  12  runs.  Runs  1-4  were  without  a  tidal  current, 
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and  Runs  5-8  had  a  moderate  steady-state  ebb  current  of 
approximately  11  cm/s  at  the  inlet  entrance.  A  strong  steady-state 
ebb  current  of  approximately  22  cm/s  at  the  entrance  was  present 
in  the  Runs  9-12. 

Wave  and  horizontal  velocity  data  were  collected  with 
electrical  capacitance  gauges  and  acoustic  doppler  velocimeters 
(ADVs),  respectively.  The  calibration  error  for  wave  height  data 
was  1  mm,  and  the  accuracy  of  velocity  data  was  within  0.5%  of 
the  measured  current.  These  laboratory  data  were  compared  to 
the  wave  model  results  to  evaluate  wave  energy  dissipation 
formulations  for  combined  depth-  and  current-induced  wave 
breaking.  Because  the  data  were  collected  only  at  handful 
locations,  the  performance  of  different  breaking  formulations  in 
WABED  was  evaluated  by  simple  statistics  like  normalized  bias 
and  root-mean-square  error. 

In  the  present  study,  the  wave  modeling  of  interest  is  at  the 
inlet  channel,  where  wave  height,  current  velocity,  and  water 
depth  were  measured  in  the  physical  model.  The  computational 
grid  covers  a  rectangular  domain  with  57  cells  in  the  cross-shore 
direction  and  12  cells  in  the  long-shore  direction.  Computational 
cell  sizes  are  30.5  x  30.5  cm2  at  the  physical  model  scale  and 


Table  3 

Model  test  conditions  for  complicated  bathymetric  coast  case 


Return  period  (year) 

1 

5 

10 

20 

100 

Significant  wave  height  (cm) 

1.86 

5.14 

5.75 

6.26 

7.26 

Significant  wave  period  (s) 

0.68 

1.00 

1.06 

1.11 

1.21 

15.3  x  15.3  m2  in  the  prototype.  Wave  gauges  were  deployed 
offshore  in  front  of  the  wave  generator,  in  the  outside  region  of  the 
inlet,  and  inside  the  inlet  area  between  two  parallel  jetties.  ATMA 
uni-directional  wave  spectrum  with  y  =  3.3  and  10  frequency  bins 
was  used  to  generate  incident  waves  at  the  seaward  boundary. 

The  water  depths  are  linearly  interpolated  from  the  isobathic 
data,  and  depths  between  the  two  jetties  are  modified  by  still- 
water  depths  at  gauges.  The  cross  section  of  the  beach  along  the 
inlet  channel  centerline  transect  is  shown  in  Fig.  1.  Current 
velocities  were  measured  at  spatially  fixed  locations.  Because  the 
detailed  two-dimensional  current  field  data  were  not  available, 
the  current  velocity  input  in  the  numerical  simulations  was 
determined  by  a  linear  interpolation  of  data.  Even  though  waves 
can  modify  the  current  field,  the  modification  of  the  current  by 
waves  was  beyond  the  scope  of  this  paper  and  therefore  it  was  not 
considered  in  the  numerical  model  simulations.  Fig.  2  shows  an 
example  of  the  interpolated  current  field  for  Run  5  and  Run  11, 
where  the  circle  and  the  square  points  denote  positions  of  wave 
gauge  and  ADV,  respectively.  Numerical  simulations  were  con¬ 
ducted  at  the  prototype  1:50  scale.  Table  1  presents  incident  wave 
parameters  and  ebb  current  conditions,  where  U3ve  is  the  averaged 
velocity  in  the  entrance  channel. 

Fig.  3  shows  the  example  of  numerical  results  for  Run  5  with 
small  peak  frequency  and  moderate  ebb  currents.  The  wave  height 
estimates  obtained  with  the  extended  Miche’s,  extended  Coda’s, 
and  Battjes  and  Janssen’s  formulas  agree  well  with  the  experi¬ 
mental  data.  In  this  example,  Chawla  and  Kirby’s  formula 
underpredicted  the  wave  height  because  the  calculated  dissipa¬ 
tion  coefficients  were  relatively  high  as  shown  in  Fig.  3(b). 


Fig.  9.  Experimental  results  of  significant  wave  heights  (normalized  by  incident 
one)  and  wave-induced  currents  for  uni-directional  waves. 


Fig.  10.  Calculated  wave  field  for  uni-directional  waves  without  currents  (using 
Battjes  and  Janssen’s  formulation). 


J.  Zheng  et  al.  /  Ocean  Engineering  35  (2008)  1090-1101 


1097 


Fig.  4  shows  the  example  of  numerical  results  for  Run  8  for  a 
large  peak  frequency  wave  with  a  moderate  opposing  current.  The 
results  using  the  extended  Miche’s  formula,  extended  Coda’s 
formula,  and  Battjes  and  Janssen’s  formula  agree  well  with  data. 
On  the  other  hand,  Chawla  and  Kirby’s  formula  substantially 
underpredicted  the  wave  height  because  the  wave  breaking 
calculated  by  this  formula  generally  occurs  further  offshore, 
producing  larger  wave  dissipation  as  shown  in  Fig.  4(b). 

Fig.  5  shows  the  example  of  numerical  results  for  Run  9  for  a 
small  peak  frequency  wave  and  strong  opposing  currents.  The 
trend  in  model  predicted  wave  heights  is  similar  to  that  for  a  small 
peak  frequency  wave  with  moderate  opposing  currents.  The  wave 
height  estimates  from  the  extended  Coda’s  formula  agree  well 
with  data.  The  calculated  wave  heights  using  the  extended 
Miche’s  formula  overestimate  the  wave  height,  owing  to  the 
relative  smaller  wave  breaking  energy  dissipation  as  shown  in  Fig. 
5(b).  The  estimates  using  the  extended  Coda’s  formula  agree  well 
with  data.  The  formulas  of  Chawla  and  Kirby  and  Battjes  and 
Janssen  underpredicted  the  wave  height. 

Fig.  6  shows  the  example  of  numerical  results  for  Run  11  for  a 
large  peak  wave  frequency  and  strong  opposing  currents.  Calcu¬ 
lated  results  show  that  the  extended  Miche’s  formula  overpredicted 
the  wave  height  while  predictions  by  the  extended  Coda’s  and  the 
Battjes  and  Janssen’s  formulas  agreed  better  with  data.  The  Chawla 
and  Kirby’s  formula  again  underestimated  the  wave  height. 

The  underprediction  of  wave  height  using  Chawla  and  Kirby’s 
formula  is  attributed  to  two  factors.  First,  the  calculated  wave  energy 
dissipation  occurs  offshore  in  relatively  deep  water  because  there  is 
no  limiting  condition  in  this  formula.  Second,  an  excessive  amount 


of  wave  energy  dissipation  occurs  because  of  unknown  value  of  y 
arbitrarily  chosen  for  an  application  that  could  be  inappropriate. 
Consequently,  it  was  necessary  to  modify  Chawla  and  Kirby’s 
formula  to  impose  a  limiting  condition  on  wave  breaking  using 
Eqs.  (9)  and  (10).  This  modification  has  been  denoted  in  Fig.  7  as 
Chawla  and  Kirby’s  formula  with  the  breaker  criterion.  A  second 
alternative  to  modify  the  Chawla  and  Kirby’s  formula  was  to  add  a 
fraction  of  the  breaking  waves  as  defined  in  Eq.  (21).  This 
modification  is  denoted  in  Fig.  7  as  Chawla  and  Kirby’s  formula 
with  a  breaker  percentage.  A  third  alternative  was  to  use  y  =  1.0  in 
Eq.  (18).  Results  of  wave  heights  calculated  from  these  modifications 
are  shown  in  Fig.  7.  By  imposing  the  limit  condition  with  Miche’s 
breaker  criterion  (the  first  modification)  gave  rise  to  a  much  stronger 
wave  convergence  (concentration  of  waves)  in  the  intermediate 
water  depth.  The  application  of  a  limit  condition  to  eliminate  excess 
dissipation  if  waves  were  not  breaking  (the  second  modification)  did 
not  substantially  improve  model  predictions  as  compared  to 
measurements.  For  the  large  peak  frequency  wave  and  moderate 
opposing  current,  using  y  =  1.0  (the  third  modification)  in  Eq.  (18) 
produced  only  slightly  better  wave  height  estimates  for  all  runs. 

Two  statistical  parameters  are  applied  to  evaluate  the  overall 
performance  of  four  different  parameterized  energy  dissipation 
formulas  in  WABED  (Table  2).  The  first  statistics  is  the  mean  value 
of  absolute  relative  error  for  the  normalized  significant  wave 
height,  defined  as 


i  M 

MAERH  =  —  V 


i  r  h,/,  ■ 

L,-( 

‘  H„3  1 

mi  1 

|  |(Wl/3)o 

(Hl/3>oJ 

1 

Hy  3 

1  mi 

(Wj/3)o 

x  100% 


(22) 


Fig.  11.  Calculated  wave  field  for  uni-directional  waves  with  currents  (using 
extended  Miche’s  formulation). 


Fig.  12.  Calculated  wave  field  for  uni-directional  waves  with  currents  (using 
extended  Goda’s  formulation). 
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where  M  =  60  (total  numbers  of  wave  height  data  available  in 
each  experiment  condition),  and  subscripts  c  and  m  denote  the 
calculated  and  measured  normalized  significant  wave  height, 
respectively.  A  small  value  of  MAERH  indicates  a  good  agreement 
between  calculated  wave  model  predictions  and  the  experimental 
data.  A  value  of  zero  implies  a  perfect  match  between  computa¬ 
tions  and  measurements.  The  second  statistics  used  here  is  the 
correlation  coefficient  (equivalent  to  the  normalized  root-mean- 
square  error)  between  calculated  and  measured  normalized  wave 
heights  (Mase  et  al.,  2005b). 

Based  on  comparisons  of  calculated  and  measured  wave 
heights  in  the  idealized  inlet  experiment  study  (Figs.  3-7  and 
Table  2),  the  extended  Coda’s  and  the  Battjes  and  Janssen's 
breaking  formulas  performed  best  for  all  wave  conditions 
investigated  herein.  The  extended  Miche’s  formula  performed 
well  for  small  frequency  waves  or  moderate  currents,  but  it 
overpredicted  wave  heights  for  large  peak  frequency  waves  with 
strong  opposing  currents.  The  Chawla  and  Kirby’s  formula 
produced  the  largest  wave  dissipation  and  therefore  under¬ 
predicted  wave  height  in  all  simulations.  For  large  peak  frequency 
waves,  this  formula  calculated  excessive  wave  dissipation  in  the 
offshore  region  for  non-breaking  waves. 

3.2.  Wave  transformation  over  complicated  bathymetry  with  strong 
nearshore  currents 

Laboratory  experiments  were  conducted  at  the  Research  and 
Development  Department  of  the  Kansai  Electric  Power  Co.  Inc., 


Fig.  13.  Calculated  wave  field  for  uni-directional  waves  with  currents  (using 
Chawla  and  Kirby’s  formulation). 


Japan  to  study  the  wave  transformation  over  a  complicated 
bathymetry  near  a  power  plant.  Fig.  8  shows  the  beach  built  in  the 
physical  model  in  a  wave  basin  of  20-m  wide  and  38-m  long,  with 
the  water  depth  contours  shown  are  in  units  of  centimeter.  The 
model-to-prototype  Froude  scale  in  these  experiments  was  1:125. 
Both  uni-directional  and  multi-directional  waves  were  generated 
by  the  wavemaker  using  the  Bretscheider-Mitsuyasu  type  of  wave 
spectrum  (Mitsuyasu,  1970).  The  wave  directional  spreading 
function  with  a  value  of  Smax  =  25  was  used  for  the  multi¬ 
directional  random  waves.  The  wavemaker  had  60  wave  paddles, 
and  each  paddle  was  30  cm  wide.  Table  3  shows  wave  conditions 
tested  in  the  physical  model. 

Waves  and  currents  were  measured  at  127  locations  in  six 
regions  with  electrical  capacitance  wave  gauges  and  electro¬ 
magnetic  current  meters  shown  as  the  small  circles  in  Fig.  8.  The 
current  velocities  were  measured  in  the  middle  of  the  local  water 
depth.  The  calibration  error  for  wave  height  data  is  within  0.1  cm 
and  the  accuracy  of  velocity  speed  is  within  1%  of  the  measured 
current  magnitude.  Both  wave  and  current  data  were  collected  for 
400  s  after  wave  paddles  started  to  move,  and  the  first  60  s  data 
were  discarded  in  the  analysis. 

For  the  case  of  100-year’s  return  period  wave,  that  is,  the 
incident  waves  at  model  scale  having  a  significant  wave  height 
and  period  of  7.3  cm  and  1.21  s,  respectively,  measurements 
indicated  strong  wave-induced  nearshore  currents  developing, 
and  the  maximum  measured  current  speed  was  up  to  25  cm/s 
for  uni-directional  waves.  These  strong  currents  caused 
waves  focusing  in  the  vicinity  of  the  concave  coasts  as  shown  in 
Fig.  9. 


Fig.  14.  Calculated  wave  field  for  uni-directional  waves  with  currents  (using 
Battjes  and  Janssen’s  formulation). 
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The  along-shore  and  cross-shore  grid  cell  sizes  were 
0.2  m  x  0.2  m  in  the  numerical  model.  The  input  wave  spectrum 
was  divided  into  10  frequency  bins  and  36  direction  bins.  The 
current  field  input  in  each  grid  cell  was  interpolated  from  the 
measured  data.  The  variation  of  water  level  due  to  wave  setup  was 
not  considered  in  these  simulations. 

Figs.  10-14  show  heights  of  the  100-year’s  return  period 
incident  uni-directional  waves  calculated  with  different  energy 
dissipation  formulas  with  or  without  the  input  current  field.  These 
results  show  that  the  wave  height  estimates  without  the  input 
current  field  are  quite  different  from  the  experimental  data.  The 
wave  convergence  in  front  of  the  concave  shoreline  is  apparent.  In 
contrast,  the  wave  height  estimates  with  the  input  current  field 
were  comparatively  better.  This  clearly  indicates  the  importance 


Measured  H1/3/  (H1/3)0 

Fig.  15.  Comparison  between  measured  and  calculated  normalized  wave  heights. 

Section  x=9.0m 

. Without  current 

-  Extended  Miche's  Formulation 

- Extended  Goda's  formulation 


y(m) 


of  the  need  to  consider  simulating  wave-current  interactions  in 
the  wave  height  estimates.  The  comparison  of  wave  estimates  for 
different  parametric  breaking  formulas  obtained  with  the  effect  of 
current  in  Chawla  and  Kirby’s  formula  was  unable  to  produce  the 
extent  or  the  intensity  of  wave  convergence  that  was  observed  in 
the  laboratory  and  prototype  settings. 

Normalized  significant  wave  height  estimates  with  different 
parametric  energy  dissipation  formulas  are  compared  to  data  in 
Fig.  15.  The  percentage  differences  between  the  calculated  and 
measured  wave  heights  are  generally  less  than  20  percent.  The 
extended  Miche’s  and  extended  Coda’s  formulas  overestimated 
the  wave  height  for  measured  normalized  wave  height  less  than 
0.75.  In  contrast,  the  Battjes  and  Janssen’s  formula  slightly 
underpredicted  the  wave  height.  Fig.  16  shows  the  example  of 
calculated  and  measured  normalized  significant  wave  heights 
along  longitudinal  and  transverse  transects  at  x  =  9.0  m  and 
y  =  6.0  m.  In  the  simulation  without  current,  the  Battjes  and 
Janssen’s  formula  was  used  to  calculate  the  dissipation.  Calculated 
results  with  the  current  effects  show  a  better  agreement  with  data 
in  the  area  of  strong  wave  convergence,  which  was  assumed  to  be 
caused  by  the  opposing  nearshore  current. 

The  mean  value  of  absolute  relative  error  for  the  normalized 
significant  wave  height  and  the  correlation  coefficient  between 
predicted  and  measured  normalized  wave  heights  are  shown  in 
Table  4.  The  sample  size  M  in  this  case  is  635.  The  two  statistics 
calculated  for  different  energy  dissipation  formulations  indicate 
that  the  mean  error  (MAERH)  in  wave  height  prediction  is  less 
than  10  percent  and  the  correlation  coefficient  is  approximately 
0.8.  The  extended  Goda’s  and  the  Battjes  and  Janssen’s  formulas 
produced  reasonably  good  estimates.  Chawla  and  Kirby’s  formula 
with  lower  limits  (y  =  0.6)  predicted  larger  wave  heights  in  the 
offshore  for  uni-directional  waves,  and  performed  rather  poorly 
for  multi-directional  waves.  For  all  four  wave  breaking  formulas, 
the  model  results  agree  better  with  data  for  incident  uni¬ 
directional  waves  than  for  multi-directional  waves.  Our  modifica¬ 
tion  of  Chawla  and  Kirby’s  formula  by  changing  y  from  0.6  to  1.0 
overestimated  wave  heights  in  the  shallow  water  where  measured 
normalized  wave  height  were  less  than  0.75.  For  the  wave 
conditions  simulated  in  this  application,  the  modified  Chawla 
and  Kirby  formula,  either  with  the  breaker  criterion  or  breaker 
percentage,  performed  better  than  the  original  formula,  producing 

Section  y=6.0m 

-  Without  current 

-  Extended  Miche's  Formulation 

- Extended  Goda's  formulation 


x  (m) 


Fig.  16.  Normalized  wave  height  comparisons  along  longitudinal  and  transverse  transects;  simulations  without  current  based  on  Battjes  and  Janssen’s  formulation: 
(a)  longitudinal  transect  and  (b)  transverse  transect. 
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Table  4 

Statistical  mean  errors  and  correlation  coefficients  for  complicated  bathymetric  coast  case 


Uni-directional  waves 

Multi-directional  waves 

MAERH  (%) 

Correlation 

coefficient 

MAERH  (%) 

Correlation 

coefficient 

Extended  Miche’s  formulation 

6.55 

0.83 

7.44 

0.79 

Extended  Goda’s  formulation 

5.74 

0.87 

5.92 

0.87 

Chawla  and  Kirby’s  formulation  (y  =  0.6) 

6.45 

0.83 

7.33 

0.81 

Chawla  and  Kirby’s  formulation  (y  =  1.0) 

6.96 

0.79 

7.68 

0.78 

Chawla  and  Kirby  ‘s  formulation  (y  =  0.6  add  breaker  criterion) 

5.79 

0.87 

5.96 

0.87 

Chawla  and  Kirby’s  formulation  (y  =  0.6  add  breaker  percentage) 

5.48 

0.88 

6.09 

0.88 

Battjes  and  Janssen’s  formulation 

5.71 

0.88 

6.12 

0.88 

similar  wave  height  estimates  as  those  obtained  with  the 
extended  Coda’s  and  Battjes  and  Janssen’s  formulas. 

4.  Concluding  remarks 

The  suitability  of  four  different  parameterized  wave  breaking 
formulas  is  examined  in  this  paper  for  implementation  in  a 
coastal  wave  transformation  model.  These  formulas  were  eval¬ 
uated  in  a  phase-averaged  spectral  model  WABED  that  incorpo¬ 
rates  the  effects  of  wave  diffraction,  reflection,  and  wave-current 
interaction.  Two  data  sets  obtained  from  laboratory  physical 
model  experiments  were  used  to  evaluate  appropriateness  of 
these  formulas  in  practical  applications  of  WABED  model.  The  first 
data  set  was  from  an  experimental  study  for  an  idealized  inlet, 
and  a  test  condition  representing  a  slack  tide  state  with  two 
steady-state  ebb  currents.  The  second  experimental  study  in¬ 
vestigated  random  wave  transformation  over  a  complex  bathy¬ 
metry,  over  which  wave  breaking  occurred,  generating  strong 
wave-induced  nearshore  currents. 

Performance  of  these  parameterized  wave  breaking  and 
dissipation  formulas  is  investigated  through  extensive  testing, 
and  results  are  provided  for  a  wide  range  of  wave  and  current 
conditions.  Statistics  of  the  predicted  wave  height  and  measured 
data  comparison  shown  in  Tables  2  and  4  revealed  that  the 
extended  Coda’s  and  the  Battjes  and  Janssen’s  breaking  formulas 
consistently  produced  reliable  and  accurate  wave  height  estimates 
as  compared  to  other  formulas.  Therefore,  these  formulas  are 
suitable  for  representation  of  wave  breaking  and  dissipation  in 
coastal  spectral  wave  transformation  numerical  models  developed 
for  coastal  inlets,  structures,  and  navigation  projects. 
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